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Abstract 

We numerically investigate the influence of intrinsic channel noise on the dy- 
namical response of delay-coupling in neuronal systems. The stochastic dy- 
namics of the spiking is modeled within a stochastic modification of the stan- 
dard Hodgkin-Huxley model wherein the delay-coupling accounts for the finite 
propagation time of an action potential along the neuronal axon. We quan- 
tify this delay-coupling of the Pyragas-type in terms of the difference between 
corresponding presynaptic and postsynaptic membrane potentials. For an el- 
ementary neuronal network consisting of two coupled neurons we detect char- 
acteristic stochastic synchronization patterns which exhibit multiple phase-flip 
bifurcations: The phase- flip bifurcations occur in form of alternate transitions 
from an in-phase spiking activity towards an anti-phase spiking activity. In- 
terestingly, these phase-flips remain robust in strong channel noise and in turn 
cause a striking stabilization of the spiking frequency. 

Keywords: synchronization, channel noise, stochastic Hodgkin-Huxley, 
delayed coupling 



1. Introduction 

Time-delayed feedback presents a common mechanism which is found in 
many biological systems including neuronal systems. Signal transmission time 
delays in neuronal systems either result from (i) chemical processes in the neu- 
ronal synapses where neurotransmitters are released and diffusively overcome 
the synaptic cleft and/or (ii) from the finite propagation speed of electrical exci- 
tations along the neuronal axon. Time delays stemming from chemical synapses 
are of the order of a few milliseconds, while the axonal conduction delays in both, 
delay-coupled neurons and autaptic feedback loops, reach values up to tens of 
milliseconds [U0]. 

As the time scale of the delayed coupling and of the neuronal dynamics 
become comparable, the delay-coupling gives raise to peculiar synchronization 
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phenomena [5]. In particular, phase synchronization phenomena in neuronal 
systems are commonly thought to be the basis for many biological relevant pro- 
cesses occurring in the brain [6l [7]. Synchrony of neurons from small brain 
regions up to large-scale networks of different cortices comes along with trans- 
mission time delays. Theoretical and computational studies on neuronal net- 
works with delay-coupling recently highlighted the occurrence of so-called phase- 
flip bifurcations [5HTU]. The ensemble activity of the coupled neurons change 
abruptly from in-phase to anti-phase oscillations or vice versa. 

With this work we research this objective by considering the influence of 
internal noise. It is an established fact that noise leads to various prominent 
effects in neuronal dynamics llj. Some typical examples that come to mind are 
stochastic resonance features [T^HTB] , and noise-assisted synchronization IT51 — 
rrS] , Within our work the intrinsic noise is due to the stochastic gating of the 
ion channels, i.e. the so-called channel noise which is inherently coupled to the 
electrical properties of the axonal cell membrane [T§H2"T] . Interestingly, it has 
been shown that intrinsic channel noise does not only lead to the generation 
of spontaneous action potentials [22], but as well affects the neuronal dynam- 
ics at different levels, namely: (i) it can boost the signal quality [THUS], (ii) 
enhance the signal transmission reliability [23] , (iii) cause frequency- and phase- 
synchronization features [24H28| and (iv) may result in a frequency stabilization 
[2T)j , to name but a few. 

The present work is organized as follows: In Sec. [2] we introduce the biophys- 
ical model. We review the standard Hodgkin-Huxley model and its generaliza- 
tions with respect to intrinsic channel noise and a delay-coupling. Numerical 
methods for simulation are introduced after that. In Sec. [3j the dynamics of a 
network of two delay-coupled Hodgkin-Huxley neurons is explored both, in the 
deterministic limit and under consideration of channel noise. As a comparison, 
we retrospect on the previous work on a single neuron subjected to a delayed 
feedback loop resulting from autapse in Sec. [4] Our conclusions are given in 
Sec. [1 

2. Biophysical model setup 

We consider a minimal building block of a neuronal network composed of two 
coupled neurons. As an archetype model for nerve excitation of the individual 
neuron, we utilize a stochastic generalization of the common Hodgkin-Huxley 
model. The stochastic generalization accounts for intrinsic membranal con- 
ductance fluctuations, i.e. channel noise, being caused by random ion channel 
gating. Moreover, we account for a delay in the coupling which accounts for a 
finite propagation time of the action potential along the axon. 

2.1. Hodgkin-Huxley-type modelling of two delay-coupled neurons 

According to Hodgkin and Huxley, the dynamics of the membrane potential 
Vi with i = 1, 2 of two coupled neuronal cells is given by [3U] 

C±Vi + G K (n) (Vi - V K ) 
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+G Na (m, h) - V Na ) + G L (Vi - V L ) = h(t) 



(1) 



Here, Vi denotes the membrane potential of the i-th cell. The stimulus I{(t) 
acting on the i-th neuron reads: 

h{t) = hMt) + IIAt) , i,j = l,2,i?j, (2) 

where the bi-directional delay- coupling of Pyragas-type |31j between the two 
neurons is assumed to be linear in the difference of the membrane potentials of 
a primary, z-th neuron at time t and a secondary, j-th neuron at an earlier time, 
t — T. The coupling thus reads: 

ir j (t) = K[V j (t-r)-V i (t)], (3) 

where k corresponds to the coupling strength and r denotes the finite delay 
time. The coupling defined in Eq. (|3| is of "electrotonic" type, i.e. we consider 
an idealized situation wherein the coupling is proportional to the difference of 
presynaptic and postsynaptic membrane potentials. This kind of coupling then 
corresponds to so-called gap-junctions which allow the bi-directional transport 
of ions and small molecules from one neuronal cell into another. Unlike the 
conductance of chemical synapses, the conductance of gap-junctions is indepen- 
dent of the presynaptic and postsynaptic membrane potential and can therefore 
be modelled by the constant coupling parameter k. Possible chemical mecha- 
nisms occurring at the synaptic cleft are assumed to be instantaneous as the 
time scale for signal propagation along the neurons axon is much larger than 
the corresponding one for the transport process in the synaptic cleft. Note, that 
the delayed stimulus in Eq. ^ results in an excitatory coupling mechanism in 
which the spiking of neuron i at an earlier time t — T time favors the initiation 
of a action potential of the other cell at time t. 

In addition to the delayed, bilinear coupling current we apply a constant 
current stimulus 2i iex t on the neurons, mimicking the common stimulus of 
the neuronal environment on the so considered two-neuron network. In ab- 
sence of the bi-directional coupling the dynamics of each neuron exhibits a 
bifurcation scenario exhibiting a subcritical Hopf bifurcation. As a conse- 
quence, the membranal dynamics displays (i) a stable fix-point, i.e. the so- 
called rest state for /j )ex t < h ~ 6.26/iA/cm , (ii) a stable spiking solution for 
h ext > I2 ~ 9.763^A/cm 2 and (iii) a bistable regime for which the stable rest 
state and a stable oscillatory spiking solution coexist, i.e. for Ii < 1^ ext < I2 
[32H3S]- In particular, for ij iex t = the membrane potential is Vrest = — 65.0mV. 

Throughout this work the membrane potentials are measured in units of mV 
and time in units of ms. For a squid giant axon, the parameters in Eq. ([I]) read 
V Na = 50 mV, V K = -77 mV, V h = -54.4 mV, and C = 1/iF/cm 2 . Further- 
more, the leakage conductance is assumed to be constant, Gl = 0.3mS/cm 2 . 
On the contrary, the sodium and potassium conductances are controlled by the 
voltage-dependent gating dynamics of single ion channels and are proportional 
to their respective numbers. In the Hodgkin-Huxley model [30], the opening of 
the potassium ion channel is governed by four identical activation gates, being 
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characterized by the opening probability n. The channel is open when all four 
gates are open. In the case of sodium channel, the dynamics is governed by a set 
of three independent and identical fast activation gates (m) and an additional 
slow, so-termed inactivation gate (h). The independence of the gates implies 
that the probability -Fk.Na of the occurrence of the conducting conformation is 
_Pk = ?i 4 for a potassium channel and P^ a = ni 3 h for a sodium channel, re- 
spectively. In a mean field description, the macroscopic potassium and sodium 
conductances then read: 

G K (n)=gr x n\ G Na (m, h) = g^ m 3 h , (4) 

where g^ ax = 36mS/cm 2 and g^ x = 120mS/cm 2 denote the maximal conduc- 
tances (when all channels are open). The two-state, opening-closing dynamics 
of the gates is governed by the voltage dependent opening and closing rates 
a x (V) and P X (V) (x = m, h, n), i.e., [30] 

n (V\ 0-0! (^ + 55) , . 

a " (y)= l-exp[-(^ + 55)/10]' (5) 

/3 n (V) =0.125 exp[-(F + 65) /80] , (6) 
_ 0.1(^ + 40) 

m[ >~ l-exp[-(V + 40)/10] ' [ ) 

P m (V)=4exp[-(V + 65)/W] , (8) 

a h (V) = 0.07 exp [- (V + 65) /20] , (9) 

^= l + exp[-(y + 35)/10] - ( 10 ) 

Hence, the dynamics of the opening probabilities for the gates read: 

x = a x (V) (1 — ar) — /3 X (V) x, x = m,h,n. (11) 

The voltage equation ([T]), Eq. Q and the rate equations of the gating dynam- 
ics Eqs. ([5])- ( [TT| ) then constitute the original, strictly deterministic Hodgkin- 
Huxley model for spiking activity of the squid giant axon. 

2.2. Modelling Channel Noise 

In this study, however, each channel defines a bistable stochastic element 
which fluctuates between its closed and open states. As a consequence, the 
number of open channels undergoes a birth-death stochastic process. Applying 
a diffusion approximation to this discrete dynamics, the resulting Fokker-Planck 
equation can be obtained from a Kramers-Moyal expansion [371 Ell ■ The cor- 
responding Langevin dynamics, interpreted here in the stochastic Ito calculus 
[3"9"j . reads: 

~x = a x [V) (1 - re) -j8 x (V) a; + &.(*), x = n,m,h. (12) 

It is driven by independent Gaussian white noise sources £ x {t) of vanishing 
mean which account for the fluctuations of the number of open gates. The 
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(multiplicative) noise strengths depend on both, the membrane voltage and the 
gating variables. Explicitly, these noise correlations assume the following form 
for a neuron consisting of N-^a sodium and -/Vk potassium ion channels: 



« m (t)f m (0) 



(1 — m)a rn + mf3. 



m 



S(t-t'), 



(13) 



N Na 
(1 - h)a h + hj3 h 

N Na 
(1 - n)a n + nj3 n 



S(t 



5(t 





(14) 



(15) 



The fluctuations of the number of open ion channels result in conductances 
fluctuations of the cell membrane eventually leading to spontaneous action po- 
tentials. These spontaneous spiking events occur even for sub-threshold, con- 
stant external current stimuli, i.e. for i^ext < I\- If the times of spike occur- 
rences are given by t n with n = 0, 1, 2, N, where N is the number of observed 
spikes, the interspike interval between two succeeding spike is T n = t n — t n —\ 
(n = 1, .., N). For the case of a single Hodgkin-Huxley neuron the distribution 
of these interspike intervals exhibits a peak-like structure with the peak located 
around the intrinsic time Ti n t r i ns i c , which is determined by the limit cycle of the 
deterministic dynamics [22] . 

The strength of the channel noise scales inversely with the number of the 
ion channels. Consequently, the threshold for excitation can be reached more 
easily with increasing the noise strength (i.e. smaller system size). In order to 
characterize the spontaneous spiking, we introduce the mean interspike interval 



With increasing number of ion channels, i.e. decreasing channel noise level, the 
mean interspike interval increases exponentially for a vanishing current stimulus 



In presence of a finite positive constant current stimulus Text, the mean inter- 
spike interval is always smaller than that for the undriven case (40] . Moreover, 
for supra-threshold driving, i.e. / cxt > I2, noise-induced skipping of spikes is 
observed. Accordingly, the channel noise docs not only favor the generation of 
spikes, but can as well suppress deterministic spiking. For intermediate constant 
current driving, i.e. I\ < I cxt < I2, for which the Hodgkin-Huxley model ex- 
hibits a bi-stability between a spiking and a non-spiking solution, channel noise 
results in transitions between these two states. 

2.3. Numerical Methods 

Our numerical results are obtained via the numerical integration of the 
stochastic dynamical system given by Eqs. (|l^-([l"2"|). Particularly, we apply the 
stochastic Euler- algorithm in order to integrate the underlying stochastic dy- 
namics 41 . An integration step At = 0.01 ms has been used in the simulations; 




(16) 



71=1 



ext — 



[22]. 
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for the generation of the Gaussian distributed random numbers, the Box-Muller 
algorithm [42: has been employed. The occurrence of a spiking event in the volt- 
age signal Vi{t) is obtained by upward-crossing of a detection threshold value 
Vq = 0. It turns out that this threshold can be varied over a wide range with 
no influence on the resulting spike train dynamics. 

To ensure that throughout all times the non-negative gating variables take 
on values solely between (all gates are closed) and 1 (all gates are open), we 
implemented numerically reflecting boundaries at and 1. Throughout this 
work we assume a constant ratio of the numbers of potassium and sodium chan- 
nels which results from constant ion channel densities. For this work we have 
assumed channel densities of 20 potassium channels and 60 sodium channels per 
fim 2 . 

In performing the numerics we initially prepare each neuron in the rest state 
voltage value V Tes t- By applying a short current pulse on one of the two neurons, 
we initialize an action potential in this neuron which later on can be echoed by 
the delay-coupling. 



3. Synchronization 

In order to investigate the temporal correlation between the spiking statistics 
of the two neurons, we apply the linearly interpolated, instantaneous time- 
dependent phase $i(f) of a stochastic spiking process of neuron i (i = 1,2) 
[THUS]; i.e., for t € [\ n ,t i>n+1 ] 

^(t) = 27rn + 27r - - k " , n = 0, 1, 2, ..JV< - 1 , (17) 

where denotes the n-th spiking of neuron i and N% is the total number of 
spike events in the dynamics of neuron i. Note, that each spike occurrence 
contributes to the overall phase with 2tt. Between two succeeding spikes the 
phase is obtained by linear interpolation. 

Note, that in accordance with this definition of the phase, the angular spiking 
frequency Ui is given by: 

Ui = lim . (18) 

t-^oo t 

This frequency and the mean spiking rate being the inverse of the mean 
interspike interval {T)i in Eq. (16) are equivalent up to a constant factor; i.e. 
we obtain, 

(19) 

Due to this equivalence we present our simulation results solely in terms of the 
mean interspike interval (T). 

In the pursuit for a measure for the synchronization of spike occurrences in 
the two coupled neurons, we consider the phase difference between the two spike 
trains, namely: 

A$(t) = ($i{t) - $ 2 (*)) mod27r . (20) 
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(a) (b) 
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Figure 1: (color online) The interspike interval (T) for the spiking of one of the two delay- 
coupled standard Hodgkin-Huxley neurons with no external constant current stimulus, i.e. 
^oxt = 0, as function of the delay time r for coupling strength re = 0.2mS/cm 2 in panel (a) 
and as function of both parameters of the delay-coupling, i.e. the coupling strength re and the 
delay time t, in panel (b). The color encoding for the phase diagram in panel (b) is shown in 
the color bar: the darker the color the shorter is the interspike interval (T). For white regions, 
the mean interspike interval (T) tends to infinity and no repitive firing is observed. 

It turns out that the phase difference AQ(t) numerically tends, after transient 
effects have faded away, for t — > oo to a value that depends on various parameters 
such as the coupling strength, the delay time, the total integration time and to 
some extent (due to inherent irregular, chaotic behavior) even on the time step 
used in our integration step. Accordingly, we observe in the long time limit that 
wi = w 2 and (T)x = (T) 2 =: (T). 

3.1. Deterministic limit 

We consider first the deterministic limit by letting A^ a — > oo and A/jc —> oo. 
In doing so, the occurrence of repetitive firing was systematically analyzed by 
varying the two coupling parameters, i.e. the coupling strength n and the delay 
time r. 

For J cxt = 0, i.e. for a subthreshold constant current driving, the resulting 
mean interspike interval (T) is depicted in Fig. [TJa). For the coupling pa- 
rameters taken within the white region of Fig. |TJb), the system relaxes to the 
non-spiking rest state. However, in the regime of the spiking dynamics a spike 
in any of the two neurons generates a subsequent spike in the other neuron. 
Consequently, repetitive, but alternating, firing can be observed for both neu- 
rons. The mean interspike interval (T) increases linearly with increasing delay 
time r, cf. Fig. [lja). In particular, 

(T) « 2(T act + r) , (21) 

where T act ~ 2ms is the activiation time which is the time between the time 
the stimulus of delayed coupling sets in and the occurrence of the stimulated 
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Figure 2: (color online) The phase difference A$, cf. Eq. l |20| , for two delay-coupled standard 
Hodgkin-Huxley neurons for 7 ox t = 0, cf. Eqs. (|T||-| |11[ |, is depicted as function of the delay 
time r in panel (a). In panel (b) the phase difference A$ is depicted within a phase diagram 
upon varying the two coupling parameters k and t. In the white region, there is no repetitive 
firing, cf. Fig. ^ In the red region the two neurons fire alternatingly, resulting in a phase 
difference of A<J> ~ it. 



spiking. The factor '2' in Eq. (21 1 is due to the alternating spiking of the two 



neurons. Note, that the time between succeeding spiking events of the network 
is the delay time plus the activation time. 

In Fig. [2] the steady-state phase difference A<E> is depicted as function of the 
coupling strength k and the delay time r. For large delay times r the neuronal 
dynamics of an individual neuron possesses after each spiking event sufficient 
time to relax back to its rest state before the delayed stimulus caused by the 
spiking event of the other neuron sets in. This results in an alternating firing 
dynamics of the two neurons. The spiking of the two neurons therefore exhibits a 
constant phase shift of ir, i.e. the spiking event of one neuron is almost perfectly 
located between two succeeding spiking events in the other neuron. However, for 
delay times that are of the order, or are even smaller than the refractory time, 
an irregular behavior of the Hodgkin-Huxley dynamics emerges, as it is to be 
expected in presence of finite delay. This in turn is reflected in the numerically 
evaluated phase difference A<I> by a noisy behavior that succeedingly smoothes 
for increasing large delay times, yielding the afore mentioned asymptotic phase 
shift of 7T. This result is corroborated with our numerics, as depicted with Fig. [2} 



3.2. Influence of channel noise 

In Fig. [3] we depict the dependence of the mean interspike interval (T) on 
the coupling parameter k and delay time r for three different levels of channel 
noise, i.e. different sets of ion channel numbers and Nk- With increasing 
noise level, i.e. decreasing number of ion channels, the sharp transition between 
the parameter regime of repetitive spiking and non-repetitive spiking is smeared 
out. 
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Figure 3: (color online) The interspike interval (T) for the spiking of one of the two delay- 
coupled standard Hodgkin-Huxley neurons with no external constant current stimulus, i.e. 
^oxt = 0, as function of the coupling parameters k and the delay time r. The mean interspike 
interval (T) is depicted for three different strengths of the channel noise: (a) strong intrinsic 
channel noise with Af[j a = 360 , TVk = 120, (b) moderate intrinsic channel noise with A^Na = 
3600 , N K = 1200 and (c) weak intrinsic channel noise with JV Na = 36000 , N K = 12000. 





Figure 4: (color online) The steady-state phase difference A$ for two delay-coupled stochastic 
Hodgkin-Huxley neurons for 7 ex t = is depicted as function of the delay coupling parameters 
k and r for three different noise strengths: (a) strong intrinsic channel noise with JVNa = 
360 , Ny^ = 120, (b) moderate intrinsic channel noise with iVjsr a = 3600 , Nr = 1200 and 
(c) weak intrinsic channel noise with Afj>j a = 36000 , = 12000. The stripes in panel (a) 

indicate noise-induced phase-flips, see text. 
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Figure 5: (color online) In panel (a) the mean interspike interval (T) is depicted as function 
of the coupling time r for constant coupling strength k, I cx t = and channel noise level 
corresponding to numbers of ion channel: A^Na = 360 and Nj^ = 120. For the same parameters, 
the phase difference A4> as function of the delay time r is depicted in (b). 



In addition, for considerable strong channel noise, distinct synchronization 
patterns emerge, indicating a frequency locking similar to the one observed for 
the autaptic case discussed in Sec. [4] below. In order to analyze the observed 
synchronization patterns in greater detail, we depict the mean interspike interval 
(T) as function of the delay time r for a fixed coupling strength k = 0.7mS/cm 2 
in Fig. [5^ a). We find that the mean interspike interval varies with the delay time 
in an almost piecewise linear manner, displaying sharp, triangle-like textures. 

The distinct peak locations of the mean interspike interval (T) can be ex- 
plained by the number of spikes that match in accordance with the intrinsic 
time ^Intrinsic a full propagation time length, given by twice the delay time. The 
mean interspike interval (T) henceforth is proportional to the ratio of twice the 
delay time and the number of spikes fitting into this very delay time interval, cf. 
Eq. (21). Accordingly, channel noise results in a stabilization of the interspike 
interval to an interval around the intrinsic time scale of the Hodgkin-Huxley 
oscillator. 



3.3. Phase-flip bifurcation 

Note, that this phenomenon complies with similar pattern appearing in the 
diagram for the phase difference, cf. Fig. [4j Furthermore, we depict the phase 
difference A$ in Fig. [5jb). These sharp transitions are accompanied by pro- 
nounced noise-induced phase-flips. At the corresponding phase-flip values at 
and 7r the spiking of the basic network changes from an in-phase towards an 
anti-phase spiking: For A$ « 0, both neurons spike simultaneously, implying 
synchronous firing, cf. Fig. [6j 
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Figure 6: (color online) The simulated time-course of the membrane potentials for the two 
delayed coupled Hodgkin-Huxley neurons is depicted: (a) for the parameters k = 0.7mS/cm 2 
and r = 8ms, the dynamics shows alternating, i.e. antiphase, spiking, (b) for k = 0.7mS/cm 2 
and r = 15ms synchronous, inphase firing of the two neurons is observed. 



4. Retrospect: Frequency locking by an autaptic feedback loop 

The observed frequency stabilization in the dynamics of a network of two de- 
layed coupled neurons shares it's origin with the frequency locking phenomenon 
in noisy neurons with an autaptic feedback-loop. When neuronal dendrites es- 
tablish an autapse, i.e. a connection to the neuron's own axon, a delayed feed- 
back loop is induced to the neuron's dynamics. Such autosynapses, described 
originally by Van der Loos and Glaser in 1972 [5] are a common phenomenon 
found in about 80% of all analyzed pyramidal cells in the cerebral neocortex of 
human brain [3]. 

Auto-synapses establish a time-delayed feedback mechanism on the cellular 
level |3j. From a mathematical point of view, autaptic connections introduce 
new time scales into the single neuron dynamics which gives raise to peculiar 
frequency looking behavior Our modelling, cf. Eqs. ([l|-([l2|, captures the 
stochastic autapse phenomena for i = j and only one neuron. 

In the limit of vanishing channel noise by letting N-Na — > oo and TVk —> oo, 
the spiking period is given by the delay time r plus the activation time T act , 
being the time needed for creating the next spike event after the delayed stimulus 
did set in, cf. Fig. [7]ja). Note, that in presence of the autaptic delay, the fixed- 
point solution of the unperturbed Hodgkin-Huxley dynamics remains stable and 
the delayed stimulus is excitatory. 

In presence of finite channel noise, however, there are two competitive mech- 
anisms at work: (i) there are spiking repetitions due to the delay-coupling and 
(ii) noise-induced generation of spikes or (iii) noise-induced skipping of spikes. 
The interplay between these mechanisms becomes evident when the distribution 
of the interspike intervals is considered. In particular, the interspike interval his- 
togram (ISIH) exhibits a bimodal structure, exhibiting two peaks, see in Fig. 6 
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Figure 7: (color online) The interspike interval (T) for the standard Hodgkin-Huxley neuron 
with an autaptic feedback loop, cf. Eqs. |T||-| |11| |. In panel (a) the dependence of the interspike 
interval (T) on the delay time r is depicted for different coupling strengths k. In the regime 
of repetitive firing the mean interspike interval (T) grows linearly with the delay time r. In 
the phase diagram shown in panel (b) the dependence of the interspike interval (T) on the 
coupling strength k and the delay time r is depicted. The color bar next to panel (b) gives 
the color encoding for the values of the interspike intervals. The white region corresponds 
to the situation for which the externally initialized spike is not echoing itself and, formally, 
(T) -> oo. 
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Figure 8: (color online) The mean interspike interval (T) for the stochastic Hodgkin-Huxley 
neuron with an autaptic feedback loop, cf. Eqs. aTn- dTSP . The channel noise level corresponds 
to iVivj a = 360 and Wk = 120. Similarly to Fig. If] (T) is shown as function of the delay time 
t in panel (a), and the corresponding phase diagram is depicted in panel (b). The color bar 
gives the color encoding of the mean interspike intervals (T). 
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in Ref. [25] . Upon specific values of the noise strength determined by the num- 
ber of ion channels and the coupling parameters k and r the bimodal structure 
shows distinct differences: Due to first mechanism the delay coupling leads to 
a significant peak around the delay time r. Via the noise-induced mechanisms, 
the channel noise results in a broad peak around the intrinsic time scale Tj ntr i ns j c . 

However, for considerable strong coupling strengths the bimodal structure 
collapses to an unimodal one and a frequency-locking phenomena takes place. 
The mean interspike interval (T) becomes bounded by a finite range of values 
around the intrinsic time scale, but still shows a striking dependence on the 
delay time r, cf. Fig [8] In particular, the mean interspike interval (T) varies 
with the delay time r in an almost piecewise linear manner, displaying sharp 
triangle-like textures, cf. Fig. |8^a) for k — 0.7mS/cm 2 and r > 10 ms. 

5. Conclusions 

With this work we have investigated the effects of intrinsic channel noise 
on the spiking dynamics of an elementary neuronal network consisting of two 
stochastic Hodgkin-Huxley neurons. In doing so, we invoked some idealistic 
simplifications such as the use of bilinear coupling with identical time delays 
and equal coupling strengths. The finite transmission time of an action potential 
traveling along the neuronal axon to the dendrites is the cause for the delay- 
coupling to the other neuron. Physically this transmission time derives from 
the finite length to the connecting dendrites and the finite propagation speed. 

A Pyragas-like delayed feedback mechanism has been employed to model 
the delayed coupling. The two basic parameters for the delay-coupling are the 
coupling strength k and the delay time r. In terms of these two parameters the 
delayed feedback mechanism results in a periodic, repetitive firing events of the 
neurons. 

Apart from this repetitive firing, the delay-coupling introduces intriguing 
time scales. Particularly, we detect a noise-induced locking of the spiking rate 
to the intrinsic frequency of the system. Consequently, the delayed feedback 
mechanism serves as a control option for adjusting the interspike intervals; this 
feature may be of importance for memory storage [I3J and stimulus locked short- 
term dynamics in neuronal systems [44] . One may therefore speculate whether 
ubiquitous intrinsic channel noise in combination with the autapse phenomenon 
is constructively harvested by nature for efficient frequency filtering. 

Moreover, the emergence of a correlation between the firing dynamics of 
the two delayed coupled neurons has been studied in terms of stochastic phase 
synchronization. The dynamics of two coupled Hodgkin-Huxley neurons ex- 
hibits noise-induced phase-flip bifurcations. At these phase-flip bifurcations the 
phase difference changes abruptly and the spiking of the neuron switches from 
an in-phase spiking to an anti-phase spiking, or vice versa. These phase-flips 
are the direct result of the influence of channel noise. The observed phase-flips 
thus may possibly assist the fact of a coexistence of various frequency rhythms 
and oscillation patterns in different parts of extended neuronal networks, as for 
example it is the case for the brain. 
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